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Abstract 

We investigate the possibility of hidden non-Abelian Local Phase symmetries in large-C/ 
doped planar Hubbard antiferromagnets, believed to simulate the physics of two-dimensional 
(magnetic) superconductors. We present a spin-charge separation ansatz, appropriate to in- 
corporate holon spin flip, which allows for such a hidden local gauge symmetry to emerge in 
the effective action. The group is of the form ST/ (2) (g> J7s(l) <8> U em (l), where SU(2) is a lo- 
cal non-Abelian group associated with the spin degrees of freedom, U em (l) is that of ordinary 
electromagnetism, associated with the electric charge of the holes, and Us(l) is a 'statistical' 
Abelian gauge group pertaining to the fractional statistics of holes on the spatial plane. In a 
certain regime of the parameters of the model, namely strong Us(l) and weak SU(2), there is 
the possibility of dynamical formation of a holon condensate. This leads to a dynamical breaking 
of SU(2) — > U(l). The resulting Abelian effective theory is closely related to an earlier model 
proposed as the continuum limit of large-spin planar doped antiferromagnets, which lead to an 
unconventional scenario for two-dimensional parity-invariant superconductivity. 
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1 Introduction 



The discovery of the quasi-planar high-T c Cuprates [0J prompted considerable theoreti- 
cal interest in two-dimensional superconductivity of magnetic origin. The strong suppres- 
sion of the isotope effect was one of the main reasons for looking for alternatives to phonon 
mechanisms. The main feature of the magnetic superconductivity was believed to be the 
fractional statistics of the excitations on the planar geometry of the materials. In two 
spatial dimensions, particles are no longer limited to Bose and Fermi statistics but can 
acquire an arbitrary interchange phase; such particles with fractional statistics are known 
as anyons. Laughlin suggested that a gas of anyons may exhibit superconductivity at 
low temperature. This idea was supported by the results of calculations in the random 
phase approximation ||, which demonstrated that a perfect gas of charged anyons with 
certain values of the statistics parameter is indeed a superconductor at zero tempera- 
ture []. This 'anyonic superconductivity' is an entirely novel phenomenon which has no 
analog in three-dimensional systems. Motivated by the role of anyonic quasi-particles in 
the Fractional Quantum Hall Effect, Laughlin went on to suggest that the charge carriers 
in the copper oxide planes of materials such as L^CuC^ and YBa2CuC>6 might also have 
fractional statistics and that superconductivity in these materials may be well described 
by the anyonic model. 

From the experimental point of view, however, there seems to be a serious drawback 
with the anyonic model as a candidate theory of high-T c superconductivity. A field the- 
oretic realisation of anyonic matter consists of fermions interacting with an abelian 'sta- 
tistical' gauge field whose dynamics is governed by a Chern-Simons term. As discussed 
in ||, this term leads to observable parity violation in an anyonic superconductor for 
which there is, as yet, no conclusive experimental evidence. In ref. || a proposal was 
made for a simple gauge-theory model which exhibits two-dimensional superconductivity 
without parity violation. In its most general form, the model consists of two species of 
massive fermions coupled with opposite signs to an abelian gauge field representing ef- 
fective spin interactions among the holon excitations. The two species have equal and 
opposite masses and hence parity is conserved overall. This theory may be shown ||, 
to arise as an approximate long-wavelength limit of an idealised model of the dynamics 
of the charge carriers in doped t — j or Hubbard models. Similar models, but only at a 
continuum theory level with no attempt to discuss the connection with semi-microscopic 
condensed-matter systems, have been proposed simultaneously in refs. j7|, |8|. 

1 For restrictions on the validity of the RPA approximation in the context of effective field theories of 
parity- violating anyonic superconductivity see ref. B . 
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The treatment in ref. || employed large-spin approximations for the antiferromagnetic 
model. This resulted in a strong suppression of intrasublattice hopping ||, which lead 
to two species of hole excitations for the bi-partite lattices used [|| |J. One eventually 
would like to argue that the same qualitative features occur for the realistic value of the 
spin, 1/2. It is the purpose of this article to attempt to formulate the above-mentioned 
effective theory and its physical consequences in a way so as to avoid the large-spin S 
assumption. 



To this end, we first review the passage from the statistical large-spin models to the 
continuum theories, and then extend the analysis to spin 1/2 models. The local phase 
symmetries that these models possess play a crucial role in this programme, and below 
we study them in some detail. What we shall show here is that the doped large-?/ Hub- 
bard models possess a local SU{2) x Us(l) phase symmetry related to spin interactions. 
This symmetry will be discovered through a spin-charge separation ansatz, which allows 
intersublattice hopping for holons, and hence spin flip. The spin charge-separation may 
be physically interpreted as implying an effective 'substructure' of the electrons due to 
the many body interactions in the medium. This sort of idea, originating from Ander- 
son's RVB theory of spinons and holons |]IU||, seems to be pursued recently by Laughlin, 



although from a (formally at least) different perspective than the one discussed here [] 1 



The effective long wavelength model is remarkably similar to a three-dimensional gauge 
model of particle physics proposed in ref. |T2] as a toy example for chiral symmetry 
breaking in QCD. In that work, it has been argued that dynamical generation of a fermion 
mass gap due to the U(l) group in 577(2) x £7(1) breaks the SU(2) group down to a 
r 3 — U(l) group, where T3 is the 2x2 Pauli matrix. From the particle-theory view point 
this is a Higgs mechanism without an elementary Higgs excitation. 



The analysis carries over to the present well, if one associates the mass gap 

to the holon condensate. The resulting effective theory of the light degrees of freedom 
is then similar to the continuum limit of describing unconventional parity-conserving 
superconductivity. 



From our point of view, the above symmetry-breaking pattern summarizes the ef- 
fects of doping on the large-?/ Hubbard model in a dynamical way. In our opinion, 
the appearance of non-Abelian gauge symmetries, as symmetries of doped antiferromag- 
nets which are broken dynamically by doping, and the analogy of the holon-condensate 
formation/superconductivity with chiral symmetry breaking in Yang-Mills theories, open 
up many possibilities for a non-perturbative (exact) treatment of such theories, including 
the role of non-perturbative effects in the superconductivity mechanism We also 
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believe that our analysis offers quantitative support to the ideas of refs. [10, 11] about 
effective 'splitting' of electrons into spinon and holons in the medium in a more general 
context. 



The structure of the article is as follows: in section 2 we discuss the doped Hubbard 
(lattice) models from the point of view of the proposed spin-charge separation ansatz and 
the associated gauge symmetry structure SU(2) x Us(l)- In section 3 we discuss the 
long-wavelength effective lattice action in the limit of strong Us(l), the dynamical mass 
generation for the holes, and the connection with (Kosterlitz-Thouless) superconductivity. 
In section 4 we present an analytical derivation of the dynamical breaking of the SU(2) — > 
£7(1) on the lattice, in the limit of strong Ug(l). Finally in section 5 we present our 
conclusions and the possible predictions of the model. 



2 Hubbard Models and Local Phase Symmetries 

2.1 Large-Spin Treatments and their Continuum Limits 

First let us briefly review the large-spin treatments of antiferromagnets || ||. In the 
absence of doping impurities, the quasi-planar materials are antiferromagnetic insulators. 
The potential importance of antiferromagnetic correlations for high-temperature (cuprate) 
superconductivity was first noted by Anderson [|10J who suggested that the correct model 



for the dynamics of electrons in the copper oxide layers was the single-band, large-?/ Hub- 
bard model. The two-dimensional Hubbard Hamiltonian is written in terms of operators, 
Cj. CT and c| CT , which annihilate and create electrons in the d x 2_ y 2 orbital at each copper 
site, 

where t is the electron-hopping matrix element, U is the strong Coulomb repulsion and 
n ia = c\ a c ia is the occupation number at each site. In the limit U — > oo, a single- 
occupancy constraint is rigidly imposed. The undoped case is described by the Hubbard 
model with half-filled band and hence the spins are the only degrees of freedom in this 
limit. To leading order in large-?/ perturbation theory, the half-filled Hubbard model is 



simply equivalent to the two-dimensional Heisenberg antiferromagnet [14 



H = J y £S i .S j (2.2) 

m 

where J = 4t 2 /U and Sj is the electron spin at site i. Thus, we see that in the infinite 
U — > oo limit J corresponds to a weak coupling. 
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The effective long-wavelength degrees of freedom of the antiferromagnet can be de- 
scribed by a 'relativistic' quantum field theory in (2+l)-dimensional spacetime. In par- 
ticular, the large- S limit of the spin-S* Heisenberg antiferromagnet is equivalent, at large 
length-scales, to the quantum nonlinear cx-model |T5| , |TBJ. The relativistic covariance of 



the effective action arises from the linear dispersion relation for long-wavelength magnons 
and the spin-wave velocity plays the role of the velocity of light in this formulation. 



Doping introduces mobile charges which hop from site to site against the antiferromag- 
netic background of the spins. The coupled dynamics of holes and spins in the doped 
system is highly non-trivial. The hopping of holes tends to disorder the spins reducing 
the antiferromagnetic correlation length and the spins also mediate interactions between 
the holes. Roughly speaking there is competition between the influence of the spins which 
favour a Neel-ordered ground-state and that of the holes which tend to form a spin liquid. 
A general conjecture is that a superconducting pairing of holes arises out of this compe- 
tition [117, This has been verified in ref. in an effective large-spin analysis. In 
that analysis a large-spin S —>■ oo has been employed, which leads to two kinds of holes, 
due to the suppression of intersublattice hopping [§, || . To be complete, below we review 
briefly the approach of ref. [H, with which we shall make contact later on. 



To this end, we first note that to describe the dynamics of holes in the model of ref. 
one implements a spin-charge separation, which is achieved by representing the electron 
operators c i)(T using a 'slave- fermion' ansatz ]B3L EDI , 



4,a = ( 2 -3) 

where ip is a Grassmann field representing the absence of a spin at a given site (hole) which 
carries the electric charge and Zj iCr is the spin degree of freedom, which can be identified 
[19] with the magnon field of the CP 1 <r-model. The ansatz ( |2.3|) carries information 



about a local gauge invariance of the model as one can perform simultaneous local phase 
rotations on ^ and z^ a fields with opposite couplings without changing c ij(T . It is this 
symmetry that is responsible for the gauge nature of the interactions between holes. The 
physical reason for such a symmetry is the restriction of having at most one electron per 
lattice site. This redundancy of degrees of freedom is the characteristic feature of gauge 
models. 



The full partition function of the model is given as a path-integral over the Grass- 
mann fields ipi and ipj as well as the CP 1 variables z, z and a M . The corresponding 
long-wavelength limit is derived by linearizing the energy spectrum about the chemical 
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potential f|, an d is given by: 



S eff = j o dr J d 2 x-\{d ll + ia tl )z\ 2 + t> a (ip-T 3 fi+~A)V a (2.4) 

where \I/ are four- component Dirac fermions, representing holes, 7 is a constant inversely 
proportional to the J Heisenberg interaction || [0| , e is the electric charge, c is the velocity 



of light in units of the hole fermi veloccity vp — 1, and A is an external electromagnetic 
field. The Dirac nature of the holes is a result of the flux-phase assumption for each 
sublattice || . The reducible four-component representation of the Dirac algebra in space- 
time is a result of doubling, and follows directly from the local sublattice structure defined 
by the antiferromagnetic order. The opposite statistical charges of the holes in different 

sublattices leads to a r 3 coupling for the gauge field, where r 3 = ^ q ^ ^ 1 ^) ' * s ^ e 

2x2 Pauli 03 matrix representation for the generator of the T3 — U(l) group || f\ 



Integrating out the electrically-neutral magnon fields, and keeping only the leading 
terms in a derivative expansion [21. 22], one obtains the low-energy effective action of the 
electrically charged degrees of freedom: 



C 



Ag' 



:f^r + V a (ip-Sr 3f 



(2-5) 



The dimensionful gauge coupling g 2 is proportional to IK 



9 



(7) 



-1 



(2.6) 



where i] is the doping concentration in the sample. In the context of the t—j model, which 
was considered in ref. ||, this coupling may be taken strong enough so as to generate 
dynamically a gap in the hole spectrum. 



The above analysis essentially postulated the existence of two holon species, by suppress- 
ing intra-sublattice hopping. This was the result of a large-spin analysis. It is the purpose 
of this article to demonstrate that qualitatively similar long-wavelength results may be ob- 
tained for spin-| doped antiferromagnets. An important tool in such an analysis is the 
study of local phase symmetries of the model, which we now turn to. We shall start with 
a review of (non-Abelian) gauge symmetries that characterize the half-filled (undoped) 
models, and then proceed to a study of the doped case upon constructing an appropriate 
spin-charge separation ansatz, extending Q2.3|) appropriately so as to allow intersublattice 

2 Here and in ref. Q we follow the terminology T3 — U(l) for the abelian group in spin or sublattice 
space generated by the 2x2 (73 Pauli matrix, so as to distinguish it from space-time groups. 
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hopping of holons. As we shall show, under the proposed ansatz, the effective Hamilto- 
nian of holon and spinon degrees of freedom is characterised by hidden non-Abelian local 
phase symmetries. However, the holon condensate breaks the non-Abelian symmetry dy- 
namically down to the abelian subgroup discussed in ref. ||, and hence one recovers 
the above-discussed Abelian model as an effective theory of the light degrees of freedom. 
Nevertheless, there are remnants of the non-Abelian symmetry structure, which manifests 
itself in the (mass) spectrum of meson-like excitations as we shall discuss in section 3. The 
presence of such excitations constitutes physically testable predictions of the spin-charge 
separation ansatz proposed in this work. 



2.2 Half-filled Spin-^ Antiferromagnets: SU(2) Gauge Symme- 
try Structure 

The large- £/ (Mott) limit of the half-filled Hubbard model with j = At 2 /U is the Heisen- 
berg model fl2.2|) . In ref. [23] it has been observed that in this limit there is a local SU(2) 



symmetry associated with the spin-1/2 algebra of the electrons. Indeed, using electron 
operators cf at a site i, corresponding to spin components up or down, a = 1,2, one may 
represent the Hamiltonian (|2.2j) as 



H= JYtf"<<*W"<&(>)i (2-7) 



with the constraint of one electron per site: 



(2.8) 



H is invariant under the usual global SU{2) transformations of the spin-| algebra, c a — > 
c fi9ai with gP an SU (2) matrix. In ref. a second SU (2) was constructed out of the 
doublet of creation operators (c3>, —c\). Combining these two doublets in a 2 x 2 matrix 

-d) (2 - 9) 

one observes that in addition to the global SU{2) transformations Xa/3 —> Xa^gJ, one 



can 



23|] define a local SU(2) by left multiplication 

Xaf3 -> hlx-,/3 (2.10) 

This local symmetry commutes with the global SU(2) mentioned above. Writing the 
global SU(2) spin operators S appearing in (|2.2| ) in terms of \ as Sex trx^X ' 7 \ with T 
denoting matrix transposition, one can easily see that the Heisenberg interaction (|2.7| ) is 
invariant under this local SU(2), which is thus the symmetry of the large-?/ Mott limit 
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of the half-filled Hubbard model. It should be stressed of course that this is not an exact 
symmetry of the Hubbard model. As shown in ref. P3]| , the very constraint ( |2.8| ) of one 
electron per site, which in terms of x variables is expressed as 

Tr X V 3 X = (2.11) 

results in a time-dependent local gauge symmetry, when combined with the kinetic term 
in the Lagrangian 

L = \T, tr x\{ij t + A 0>i ) X i - H (2.12) 
where Aq^ acts as a Lagrange multiplier implementing the constraint, and it may be 



thought of as the third (temporal) component of the local SU{2) gauge field P3| . Such 
gauge symmetries appear as a general property of the Gatzwyler projection of one electron 
per site, due to the fact that such projections are associated with a sort of particle-number 
conservation. This local gauge symmetry connects various mean field limits of the half- 
filled Hubbard model [E3 . 



To understand the formal meaning of the above symmetry, we return to the CP 1 a- 
model, which is supposed to describe the low-energy physics of the half-filled Hubbard 
model in a bosonized framework for the spin excitations. We recall that upon resolving 
the constraint zz = 1, with z = (zi,z 2 ) a complex SU(2) doublet with boson statistics, 
the z field can be written as a 2 x 2 unitary matrix: 



zi - z 2 

Z2 Z\ 



exp(i£ a a a ) (2.13) 



were a a ,a = 1, 2, 3 are the Pauli generators of SU (2), and the real fields £j are dynamical. 
The gauged cx-model action in this representation reads 

S z = J d?x%Hr\(d^ - tgB^)z\ 2 (2.14) 

where 70 is a bare coupling constant. In this representation one is free to gauge the full 
SU(2) local gauge group in the cr-model action, in which case = B^cr a , a = 1, 2, 3, or 
its Abelian U(l) subgroup B^ = B 3 cr 3 . The action ( |2.14| ) reads 



S= I d 3 x^[ Y,m«) 2 + 9 2 Bl + Y,B a ,(-2gd,Q ] (2.15) 

J a a 

Technically the above representation separates the Goldstone modes from the rest of the 
fields relevant at low momenta |2S|. The resolution of the constraint implicit in ( |2.13| ) 
results in a standard mass term for the gauge field B, instead of the quartic coupling 
zB 2 z. 
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The possibility of gauging the full SU (2) group in the a model is equivalent to the local 
SU(2) symmetry of the Heisenberg action (|2.2|) found in ref. p3|, given that at half-filling 



only spin excitations (magnons) exist fL0 |. Of course, the equivalence is understood in 
terms of bosonization, which in 2 + 1 dimensions, unlike 1 + 1 dimensions, cannot be 
expressed in a closed form, but only as an effective derivative expansion. 

2.3 Doped spin-| Antiferromagnet and Non-Abelian Gauge Sym- 
metry Structure 

Doping is expected to break the SU(2) symmetry between creation and anihilation pairs 
of electron operators. Naively speaking, a spatial hopping term of the form c* a jC a j does 
not seem to be invariant under the local SU(2) (|2.10| ). Away from half-filling one would 
expect that only a local U(l) can survive, which in view of our spin-charge separation 
ansatz ( |2.3|) seems to be the Abelian subgroup of SU(2) associated with r 3 . This local 
subgroup is the one gauged in the CP 1 a model, and also the one associated with the 
(Berry phase) term describing static holes in the model of ref. . In this article we shall 
present a dynamical scenario by which the above symmetry breaking is achieved. The 
scenario will be remarkably similar to a three-dimensional particle-physics toy model for 
chiral symmetry breaking in QCD ||12 |. 



The key point is to try to uncover the local SU{2) symmetry in the doped case by 
generalizing the spin-charge separation ansatz ( |2.3| ). We seek a representation of the spin- 
charge separation that will allow spin flip, but would still treat the holons as 'blind' to 
the electronic sublattice structure. To this end, we propose to represent the holon degrees 
of freedom as two-component spinors in a two-dimensional 'colour' space, representing 
Dirac spin components, (ipi,ip2), whilst the spin excitations are represented by the CP 1 
doublets (zi,Z2) living in the same 'colour' space. However we amend our construction 
with a spin-flip operation, which, for the z-magnon degrees of freedom is represented by 
the conjugate doublet (—Z2,~zi)- Thus the electron anihilation operators can be expressed 
as 

<* = fa ^){ Z Z \) c * = fa ^(^f) (2 ' 16) 

while the correspondiing creation operators can be obtained by -c\, c\, with t denoting 
hermitean conjugation. We believe that this ansatz captures the qualitative features 
behind the RVB idea of Anderson |]10| on spinon and holons. Essentially ( |2.16| ) implies 



that to anihilate an electron with, say, spin up one has to remove all the components of 
the spin. The spin-charge separation ansatz implies that to some extent the holes should 
be 'blind' to the spin of the electron (sublattice structure of the antiferromagnet). This 
is correctly captured in (|2.16|) , since the hole-'spinors' in colour space are the same for 
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both electron components, whilst the (magnon) ^-doublets differ by a spin-flip operation 
defined above. 



Technically, it is convenient to combine the creation and anihilation operators, following 
the treatment of the half- filled case ( |2.9|) . To this end, we propose that for the large-U 
limit of the doped Hubbard model the following spin-charge separation ansatz occurs at 
each site i: 

i - ( Cl ° 2 \ (fa fa \ ( %i -z 2 \ , 

where the fields z a ^ obey canonical bosonic commutation relations, and are associated with 
the spin degrees of freedom, whilst the fields ipa,i, a = 1,2 have fermionic statistics, and 
are assumed to create holes at the site i with spin index a. They obey the anticommutation 
relations: 

{lpi, a , ^j,fs} = SijSaP {lpi,a, 1pj,p} = {^la, V^} = (2.18) 



The ansatz ( |2.17| ) has spin-electric-charge separation, since only the fields ip carry elec- 



tric charge. From now on, we shall refer to ip a as the 'holons', and to ~z a as (bosonized) 
'spinons'. The ansatz ( |2.17| ) is an obvious generalization of ( |2.3| ) if one allows intersub- 
lattice hopping. 



It worths noticing that the anticommutation relations for the electron fields c a ,ct, do not 
quite follow from the ansatz Q2.171 ). Indeed, assuming the canonical (anti) commutation 
relations for z (ip) fields, one obtains from ( |2.17| ) 



{ci,i,C 2 ,j} ~ 2^1,^2,1*1 

{ c \,i,cl,j} ~ wIMAj 

{ c l,i>4j} ~ { c 2,ii c l,j} ~ 

{ c a,i, c a ,j} ~ $ij [ z h^i,P + VfeVfe], « = 1, 2 no sum over i, j (2.19) 

(9=1,2 

To ensure canonical commutation relations for the c operators therefore we must impose 
at each lattice site the (slave-fermion) constraints 

^1,^2,1 = fa 2 ,ifal,i = > 

J2 [zi,^i,p + MU = 1 ( 2 - 2 °) 

(9=1,2 

Such relations are understood to be satisfied when the holon and spinon operators act on 
physical states. Both of these relations are valid in the large-?/ limit of the Hubbard model 
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and encode the non-trivial physics of constraints behind the spin-charge separation ansatz 
( |2.17| ). They express the constraint at most one electron or hole per site, which charac- 



terizes the large-?/ Hubbard models we are considering here. From the above analysis, 
therefore, it becomes clear that the ansatz fl2.17| ) does not characterize a generic Hubbard 



system, but only the appropriate large-?/ limit, where the constraint of one electron per 
site is valid. As we shall discuss in section 4, both of the above constraints ( |2.20| ) are 
consistent with the mass spectrum of the effective long-wavelength theory obtained from 
dynamical generation of a fermion condensate. 



Now let us look at the symmetry structure of the spin-separation ansatz ( 2.171) , which 



in view of the previous analysis coincides with the symmetry structure of the effective 
large-?/ Hubbard action. First, it appears to have a trivial local SU(2) symmetry, if one 
defines the transformation properties of the z fields to be given by left multiplication with 
the SU(2) matrices, and those of the ijy a * matrices by the left multiplication (|2.10|) . In 
this representation, the gauge group SU (2) is generated by the 2x2 Pauli matrices. 



The ansatz ( 2.17 ) possesses an additionallocal Us(l) 'statistical' phase symmetry, which 
allows fractional statistics of the spin and charge excitations. This is an exclusive feature 
of the three dimensional geometry. This is similar in spirit, although implemented in an 
admittedly less rigorous way, to the bosonization technique of the spin-charge separation 



ansatz of ref. [26], and allows the alternative possibility of representing the holes as slave 



bosons and the spin excitations as fermions. 

In addition, as a consequence of the fact that the fermions ip carry electric charge, one 
has an extra U em (l) symmetry for the problem. 

To recapitulate, the above analysis, based on the spin-charge separation ansatz ( |2.17| ) 
which allows spin flip, leads to the following local-phase (gauge) group structure for the 
doped large-?/ Hubbard model: 

G = SU(2) x U s (l) x U em (l) (2.21) 

where the second U em (l) factor refers to electromagnetic symmetry due to the electric 
charge of the holes. This symmetry appears as a hidden symmetry of the effective holon 
and spinon degrees of freedom obeying the ansatz ( |2.17| ). 



The presence of the Us(l) 'statistics' changing group factor will be crucial in our anal- 
ysis. As we shall discuss in the next section, in its strong coupling limit it can generate a 
mass gap [E7|, 128], EU[ for the fermionic holon fields ip, which for each hole component breaks 
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parity, thereby producing a statistics changing dynamical Chern-Simons term. However, 
due to the even number of fermionic species there is no overall parity violation in the 
model |J. Note that, since this statistical gauge field couples also to the z fields, their 
statistics will be affected as well. 



2.4 Effective Hamiltonian of the doped Hubbard antiferromag- 
net 

Next, we focus our attention in showing that the various terms in the action be express- 
ible in terms of the Xa/3 variables, which would imply that the symmetries of the large-?/ 
doped Hubbard model action, are the symmetries of the ansatz fl2.3|) . 



To this end, we first study the hopping term of the dopped hamiltonian, which broke 
explicitly the local SU(2) symmetry (|2.10 ) of the electron operators c a , cL. Let us rewrite 
this term in terms of Xa/3 variables: 

Hhop = ~ tii C ]y.A C a,fi = ~~ tij {Xi.a-yXj^a + Xi,ay( (J 3)-y/3Xj,f3a] (2.22) 
(ij) (ij) 

where 03 is a 2 x 2 Pauli matrix, and summation over the spin indices is implied. In terms 
of the spin and charge excitations, appearing in ( |2.17|) , then, the hopping term may be 
written as 



{ij) 



and is trivially \ocal-SU (2) symmetric. 



(2.23) 



To complete the analysis we should also look at the interaction terms. The Heseinberg 



(2.24) 



term (|2.7| ) can be written in the following convenient form [23 

h = -\j tr b(ix]xjxl) 



which can be linearized in terms of the fermion bilinears if one introduces in the path 
integral a Hubbard Stratonovich field A^, in a standard fashion. The result of the lin- 
earization is: 

H = £ tr[(8/J) + (xlA ijXj + h.c.)} (2.25) 

<ij> 
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We then employ the ansatz Q2.17 ), and perform a Hartree-Fock (mean field) approximation 
for the bilinears: 



< ZiZj >= \At\VijUij ; 

< ^(-Ujil + <r 3 ) + Aij)ipj >= \A % \VijU k 



;2.26) 



where, according to the previous discussion, we have used the fact that the link variables 
are SU{2) x Us(l) group elements, due to the specific transformation properties of the 
variables z and ip. In the above notation V is the SU(2) part and U denotes the Abelian 
Us(l) group element. The amplitudes \Ai\, i = 1,2, of the link variables are assumed 
frozen, as usual. By an appropriate normalization of the z and ip fields, this amplitude is 
common for both link variables, 



K 



(2.27) 



According to the discussion of ref. the amplitude K is proportional to the Heisenberg 
exchange interaction J = 4t 2 /U, with t the hopping parameter, and also to the doping 
concentration in the sample |19|. We shall come to this issue later on. 



The result of the Hartree-Fock approximation, then, for the combined hopping and 
interaction terms in the Hamiltonian is: 



H H F = E tr [( 8 M) A lA* + + °3) + Ay)(^' < ZjZi > 4) 



+ 



<ij> 



E tr [ z i < ^li-tiji 1 + °3) + > Zj + h.c. 



(2.28) 



and using ( |2.26D , (|2.27|) one opbtains: 

<ij> 

E tr [KziVijUijZj] + h.c. 

<ij> 



+ 



(2.29) 



This is the effective field theory lattice action we propose to describe the dynamicsl of the 
large U Hubbard model. It is understood the the constraints ( j2.20| ) should be taken into 
account, to complement the description. It is important to note that the 'fermion' fields ip 
are 2x2 matrices in the above representation. Notice also that the term ty (1 + cr 3 ) + Ajj 
transforms covariantly under a global U(l) symmetry generated by the Pauli matrix cr 3 . 
This global U(l) symmetry acts on the electron operators x% as Xi ~^ Ux%, with U = e ld , 
9 a global phase. The z-dependent (magnon) terms yield, in the continuum, the CP 1 
a- mo del lagrangian ( |2.14| ) ]T§| • 
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In the large-C/ Hubbard limit, we are considering here, one has the folowing order of 
magnitude estimates: 



J = At 2 /U t ~ U Vmax ■ Vmax « 1 (2.30) 

where r] max is the maximum doping concentration of the sample, above which supercon- 
ductivity is destroyed. For underdoped cuprates one may consider the case r\ max « 1. 
In this limit, one observes from (|2.29|) that the Gaussian fluctuations of the variable A^- 
are of order 0(J/\tij\), and hence suppressed, as compared to the hopping term t^. This 
means that one may approximate £;.,•( 1 + cr 3 ) + Ajj ~ £^-(1 + er 3 ). Considering the usual 
case with = t for every one may absorb such terms into an appropriate rescaling 
of the fermion fields ip. This will be understood in what follows. However, we stress once 
again that in the case of finite-?/ Hubbard models, one should consider the effects of the 
Gaussian variable A^ in the lagrangian ( |2.29| ). This will be left for future work. 



The conventional lattice gauge theory form of the action is derived upon integrating 
out the magnon fields, z, in the path integral. As discussed in |22, 21] , the result of such 
a path integration of the magnon fluctuations around the mean field yields appropriate 
Maxwell kinetic terms for the link variable VijUij, which are the dominant terms in a 
low-energy derivative expansion. The constraint of at most one electron per lattice site 
in (|2.20| ) is crucial in such a derivation, since its implementation through a Lagrange 
multiplier field, a, results in a 'mass' term for the magnon fields z, in the way explained 
in ref. [pl| . The effective Maxwell terms in the continuum are of the generic form: 



S kin oc / d 3 x— !=(F 



pi/ Q UV I 



(2.31) 



where F uu , Q nv denote the Us{l) and SU(2) field strengths respectively, and <To is a 



vacuum expectation value of the Lagrange multiplier field a. 
renormalization-group analysis yields 



An elementary one-loop 



a = M - 4nK_ 



;2.32) 



where M is a transmutation mass, and Kr is the 'renormalized' K coefficient of the CP 1 



part of the action ( |2.29|) . From the analysis of ref. |19[ we may infer that Kr oc J77, with 
J the Heisenberg interaction, and 77 the doping concentration, which for lightly-doped 
cuprates is 77 << 1. This implies that the order of magnitude of the coefficient of the 
Maxwell term Q2.31|) , resulting from the z- integration in a derivative expansion, is set by 
the Heisenberg exchange interaction field strength J. The conventional three-dimensional 
gauge coupling g 2 , of dimensions [mass], is related to K R by the simple relation (|2[ 



1/g oc K R oc J "77 



;2.33) 
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Thus, from ( p. 30 ) one obtains for the dimensionless coupling g 2 a, with a the lattice spacing 
of the antiferromagnetic Hubbard model: 



The magnitude of this coupling depends on the way the limit U — > oo is taken. Taking 
the limit of £/ — > oo such that Uarjri^ nax » 1 one obtains a smaZl /?, i.e. strong coupling 
for the Us(l) group. The limit of small (3 is crucial for the symmetry breaking patterns 
of the non-abelian SU(2) group, as we shall discuss in section 3. 



We now remark that on the lattice the kinetic (Maxwell) terms ( |2.31 ) are given by 
appropriate plaquette terms of the form: 

E [Psu(2)(l - TrV p ) + (3 Us(1) (l - TrU p )} (2.35) 

v 

where p denotes sum over plaquettes of the lattice, and Pu s (i) = A, Psu(2) = ffa = 4/3i 
are the inverse square couplings of the Us(l) and 577(2) groups, respectively. The specific 
relation between the SU(2) and Us(l) couplings is due to the appropriate normalization 
of the generators of the groups. 



At this point it is worthy of remarking that for certain Schwinger- Dyson treatments of 
dynamical symmetry breaking || Q a large- N treatment is desirable, in which case one 
assumes that the spin SU(2) group is replaced by SU(N) with N large enough. In that 
case the non-Abelian coupling is related to the Abelian one through 

Psu(N) = 2N(3 Us{1) = 2N(3 1 (2.36) 

This implies that, even in the case of strong Us(l) coupling, j3\ — > 0, the large-A^ (large 
spin) limit may be implemented in such a way so that (3su{N) is finite. This is the limit 
of the analysis of ref. [Q. We shall discuss this case in section 3, where we shall make 



contact with the results of ref. 0, where such a large- N treatment had been assumed. 

Above we did not write explicitly the chemical potential term \i Y^i, a c l a °ia which de- 
termines the doping concentration in the sample. This term is also expressed in terms 
of the x variables, and essentially has the form of Q2.22 ) but for i = j, which again 



may be expressed in a gauge invariant way upon using the ansatz ( p. 17 ). In deriving 



long-wavelength continuum limits, one linearizes the energy spectrum about the chemical 
potential 0, [| . For most of our discussion below we shall not write explicitly such terms, 
as they do not affect the symmetry structure of the theory. 
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2.5 Spinor Structure for Holons and Symmetry Breaking Pat- 
terns 



Before closing this section we would like to remark that, as a result of the 2x2 matrix 
structure of the fermion fields if) in ( |2.29| ), one may actually change representation of the 



SU{2) group, and, instead of working with 2x2 matrices, one may use a representation 
in which the fermionic matrices if} a p are represented as four- component vectors (in 'colour' 
(spin) space): 

Vctf->* + = (Vi ~4 ^2 4i) (2-37) 



It is easy to see that in this representation the SU (2) group is generated by the following 
matrices: 

n = 73=(° J).^=7b = <(_° 1 1 ), 7* = A = rySTO = ( ^ ) 

(2.38) 

where the substructures are 2x2 matrices. This is the SU(2) representation used in ref. 
[ |T2"[ in the context of three-dimensional toy models for chiral symmetry breaking. Remark- 
ably, the same type of symmetry arises in our context between creation and annihilation 
operators of holon pairs in the spin-charge separation ansatz ( ^.17] ). 



In the analysis of ref. |12[ , to be discussed in the context of the present model in the next 
section, the statistical group Us(l) group is responsible for the dynamical generation of a 
parity conserving mass < >. In terms of the dynamical variables describing creation 
and annihilation of holons, if), if>^ respectively, the parity conserving mass depends on the 
holon condensate. To see this, it is convenient to split the four-component spinors ( |2.37| ) 
into two-component ones: 

*1= (Vi ~4), *t= (>2 4) (2.39) 



In this representation the two-component spinors ^ fl2~39D will act as Dirac spinors, and 



the 7-matrix (space-time) structure will be spanned by the irreducible 2x2 representation. 

The Dirac conjugate field \1/ may be identified directly with the hermitean conjugate 
fields W in terms of holon operators. This is due to the fact that in a path integral 
over the holon fields, the conjugate fields ip^ can be considered as independent degrees of 
freedom || H> HU- I n this representation, the local SU{2) gauge group is generated by 



the familiar 2x2 Pauli matrices a a , a = 1,2,3. The parity transformation is defined as 
^! — > 0-1^2, ^2 ~~ > Ci^ij which in terms of the (microscopic) holon operators ipi, i = 1,2, 
reads: ifj\ — > 4^2 ~^ ~4i ■ With these in mind, it is straightforward to observe that 
the parity-conserving mass term \I/ 1 \I/ 1 — ^2^2 can be related to the holon condensation 



< - ^ 2 ^2 >= -2(^iVi - 4M (2.40) 
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where we took proper account of the anticommutation relations ( 2.18Q among the grass- 



mann tpi >a , a = 1,2. The terms < ip^ipa >,a = 1, 2, are holon condensates. Notice, in the 

same context, that the parity violating mass term < \^ 1 \^ 1 + \& 2 ^2 > equals an irrelevant 
constant, which may be subtracted. This result is consitent with the generic energetics 
arguments that disfavor dynamical generation of a parity-violating mass in vector like 
theories with even flavour number |3(| . 



The formation of holon coondesate due to a statistics changing 27s(l) group is similar in 
spirit to the approach of ref. M in the context of the anyonic superconductivity. However, 
as mentioned above, in our case, due to the four-component structure of the fermions, 
there is an even number of fermionic species and hence no overall parity violation. More- 
over this mass gap is not a singlet under SU (2), as we shall discuss in the next subsection, 
but transforms as a triplet |I2"|| , thereby breaking SU(2) down to its r 3 — 2/(1) subgroup. 



This is the r 3 — 27(1) symmetry of the alsatz ( |2.3| ), leading to the effective action ( |2.5|) . 
This provides a sort of dynamical breaking of the local spin SU(2) group as the result of 
introducing holes into the system. 

The breaking of the 527(2) symmetry down to its Abelian r 3 subgroup admits the 
(physical) interpretation of restricting the holon hopping effectively to a single sublattice. 
In a low-energy effective theory of the massless degrees of freedom this reproduces the 
results of ref. || |9| . This scenario can be readily seen by using the four-component spinor 
representation (|2.37|) . Clearly the two off-diagonal generators of the 527(2) group Q2.38 



73 and 75, corresponding to the gauge bosons acquiring masses dynamically due to the 
holon condensate, mix the two sublattices in the notation of ref. ||, ||. Indeed, from (|2.39|) 
it follows immediately that if a holon of spin, say, 1 is created at a site i, these generators 
would connect it to the destruction of a hole with spin 2 in the neighboring sublattice. On 
the other hand, the generator A of the unbroken r 3 — 27(1), is block diagonal, thereby not 
mixing the sublattices. The intrasublattice hopping in this approach is then suppressed 
by the mass of the gauge bosons. We are considering here the limit of infinitely strong 
2/5(1) jr^J. In such a limit the intra-sublattice hopping is completely suppressed, since 
the mass (which is proportional to the infinite condensate) is infinite |TjJ . This situation, 



therefore, describes static holes. Hole hopping is allowed for strong but finite couplings, 
in which case the holon condensates and masses are finite. 



We shall devote more discussion on the phase diagram of the theory, and its comparison 
to that of ref. |J, in the next section. We would like to close this section by noting 
that, in the context of microscopic models of the form fl2.29Q , dynamical formation of 
holon condensates, and hence destruction of antiferromagnetic order, would occur above 



a critical doping concentration [31]. To quantify the above results on symmetry breaking, 
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therefore, one needs proper lattice simulations of these models. This is left for the future. 



3 Long-wavelength limit of the spin-| doped antifer- 
romagnet 

3.1 Derivation of the Long- Wavelength Hamiltonian 



We now proceed in the long wavelength limit of (|2.29 ), (|2.35|) , in the spinor representa- 



tion for the holon fields, discussed in subsection 2.5. To this end, we assume -following 
the analysis of ref. 0- a non-trivial flux-phase for the gauge field £/s(l). This is crucial 
in yielding a Dirac form for the hole effective action [^, 23, |^]. The long- wavelength con- 



tinuum limit is then obtained in a similar way as in the abelian case of ref. || [L9|], at low 
energies, by linearizing about a specific point on the fermi surface []. Due to the assumed 
flux-phase- 7r background for the gauge field Us(l) one gets for the hopping (kinetic) terms 
of the two-spinors ( |2.39|) (ignoring interactions for brevity) : 

L U n ~ E(- 1 ) ro+ """ +rM " 1 ^W* C ( r + A) + h.c. (3.41) 

where c = 1, 2 is the colour index, not to be confused with the space-time (Dirac) index. 
The factors (— l) r °+- yield a phase e 47r = (—1) per lattice plaquette, and this result is 



produced in our case by the Us (I) flux phase background M. As discussed in ref. [32 



the form ( |3.41| ) corresponds to a Dirac form for the kinetic terms of the fermions ^ upon 



making an (inverse) Kawamoto-Smit transformation [33]: 



* c (r) = 7o ro • • -7^cM * c (r) = * c (r)( 7 lp ■ ■ ■ (llY° (3-42) 

where \l/ are two-component Dirac spinors, carrying 'colour'. We stress once again that 
the colour structure is up and above any space-time (Dirac) structure. Notice that in 
such a picture fermion bilinears of the form ^ i c ^> i c i (i=Lattice index), for instance the 

condensate ( |2.40| ), are just due to the Clifford algebra {7^ , 7^} = — 25^ u and 

(anti)hermiticity properties of the 2x27 matrices on the Euclidean lattice. This is useful 
to have in mind when we study the spectrum of meson states in section 4. 



3 In what follows we shall ignore, for simplicity, the shape of the fermi surface || and therefore deal 
with conventional relativistic lattice models. Of course, this will not be the case in a realistic condensed 
matter system, where there are known to be large fermi surfaces for holes. The relativistic nature may be 
accurate in superconducting models with nodes on their fermi surface, when linearization about a node is 
performed. However, for our purposes in this work, which are a study of the generic symmetry-breaking 



patterns of the local group (2.21), their physical implications for superconductivity, and the connection 



with the results of ref. Ifll, such relativistic models will be sufficient. 
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In what follows we shall make use of the above-mentioned (irreducible) 2x2 represen- 
tation in both the colour and space-time indices on the lattice. According to the above 
discussion, then, upon ignoring for the moment the electromagnetic interactions of holes, 
one obtains the following effective low-energy lattice action for the holon fields, originating 
from (pD,dH), 

+P 1 £(1 - trU p ) + fa XXl - trV p ) (3.43) 
v v 

where \i = 0,1,2, = exp(i0i !fl ) represents the statistical Us(l) gauge field, Vi jfl = 
exp{ia a B a ) is the SU{2) gauge field, and the plaquette terms are obtained, at low en- 
ergies, as a result of the z-magnon integration p2[ [21[] Q. The fermions \I/ are taken to 
be two-component spinors, in both Dirac and colour spaces. The quantity K, is propor- 
tional to the holon hopping matrix element, which in turn depends [[H], |6|] on the doping 
concentration, as stated eralier Q2.33 ). According to the discussion following ( 2.36 ) 



m a 



large-spin (SU(N), — > oo)) treatment, as in ref. 0, the coupling constant j3 2 — > oo, and 
hence the non-Abelian-gauge-group sector of the model is weakly coupled in this limit. On 
the other hand, the coupling of the statistical Us(l) is considered as strongly coupled, in 



the limit U — > oo. It is known, from either lattice results [28], or semi-analytic Schwinger 



Dyson (SD) type of analyses |£7, |2£|, that dynamical mass generation in a U(l) theory in 
three space-time dimensions occurs only for strong coupling, i.e. for values of the gauge 
coupling that are larger than a given critical value. This mass will break the SU(2) gauge 
group dynamically. This will be discussed in detail in section 4. 



The above limit has been studied in ref. ]12| , where the model Q3.43| ) has been used 
as a toy model for studying chiral symmetry breaking patterns of QCD. Remarkably, 
as we have described above, this model can also be used to describe the physics of the 
spin-charge separation of strongly correlated electrons in a doped Hubbard model in its 
large-t/ limit. In this analogy the holon fields ipa/3 behave like the 'quarks' of QCD, which 
are thus viewed as substructures of the physical electron Xap- It seems to us that this 
point of view is similar in spirit to that pursued in the context of anyonic models by 



Laughlin ||TT|. However, we should stress that from our point of view this 'splitting' is 
viewed as a many-body effect for the holon dynamics in such systems, and hence we do 
not ascribe to it any further significance. 



4 



We would like to mention that, technically, in order to study dynamical formation of fcrmion con- 



densates on the lattice using Monte-Carlo studies as in ref. Jl2[ |, one should add to the action (3.43) a 
bare mass term mo 2, ^ta^i, and take the limit too — > only at the very end of the computations. 
This will be irrelevant for our purposes here. 
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3.2 Symmetry Structure in the Continuum 



It will be instructive to study first the symmetry structure of the model ( |3.43| ) in the 
continuum, following the analysis of ref . . This will help the reader understand better 



the interplay between the irreducible (2 x 2) and the reducible (4 x 4) representations of 
the Dirac and colour (gauged-chiral symmetry) groups. To this end, we first note that 
the continuum limit of the model ( |3.43| ) is described by the lagrangian JT2] : 

C = -\{F^f - \{Q^f + - (3.44) 



with = — ig±a^ — ig20~ a B a ^ and F^^Q^ are the corresponding field strengths for 
the abelian (statistical) gauge field and the spin SU(2) gauge field B" The parity 
conserving bare mass m term has been added by hand, as mentioned above, to facilitate 
Monte-Carlo studies of dynamically generated fermion masses as a result of the formation 
of fermion condensates < > by the strong Us (I) coupling. The m = limit should 
be taken at the end. 



To understand better the nature of this SU(2) gauge symmetry, it is instructive to look 



first at the global SU (2) group, whose gauging produces the action (|3.44j) . To this end we 
observe that the 7 M , /i = 0, 1,2, matrices, which span the reducible 4x4 representation 
of the Dirac algebra in three dimensions in a fermionic theory with an even number of 



fermion flavours, assume the form [27] : 



7 



7 




? 




(3.45) 



where a are 2x2 Pauli matrices and the (continuum) space-time is taken to have 
Minkowskian signature. As well known |27]] there exists two 4x4 matrices which an- 
ticommute with 7^,/i = 0, 1, 2: 



73 = ^ 1 ) ' 75 

where the substructures are 2x2 matrices, 
symmetry for the massless-fermion theory 











(3.46) 



These are the generators of the 'crural' 



(3.47) 
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Note that these transformations do not exist in the fundamental two-component repre- 
sentation of the three-dimensional Dirac algebra, and therefore tha above symmetry is 
valid for theories with even fermion flavours only. 



The set of generators {1,73,75,27375 = A} form [Tj a global SU{2) x U(l) symmetry. 
The identity matrix 1 generates the U(l) subgroup, while the other three form the SU(2) 
part of the group. The currents corresponding to the above transformations are [|l^] 

Jl = ^ 7m ™ T = 73, 75, i7s7s (3.48) 

and are conserved in the absence of a fermionic mass term. It can be readily verified that 



the corresponding charges Qr = J d 2 x^T^/ lead to an SU (2) algebra [[12 



[Qs, Q 5 ] = 2iQ A [Qb, Qa] = 2iQ 3 

[Qa,Q 3 ] =2*Q 5 (3.49) 

If a mass term is present then there is an anomaly 

d^Jl = 2m^r^ (3.50) 

while the current corresponding to the generator 1 is always conserved, even in the 
presence of a fermion mass (TT 

The bilinears 

Ai = ^73^, A 2 = "*7 5 ^, A 3 = 

B lfl = ¥>7s*, B 2il = ¥>75*, = * 7m A#, /x = 0, 1, 2 (3.51) 

transform as triplets under SU(2). The SU(2) singlets are 

A A = ^A^, B^ = (3.52) 

i.e. the singlets are the parity violating mass term, and the four-component fermion 
number. 



In two-component notation for the spinors the above bilinears read W2 



Al = -^1*2-^2*1], A 2 = *l* 2 + *2*l, ^3 = ^1^1 - *2^2, 

Ai = #i*i + #2*2, £ 4 , M = (7^1 + ¥2^*2 , /I = 0,1, 2 (3.53) 
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with denoting two-component Dirac spinors. For later convenience we have passed 
onto a three-dimensional Euclidean lattice formalism, in which \l/ is identified with 
c.f. (|2.39|) . In this convention the bilinears (|3.53| ) are hermitean quantities. It is this 
Euclidean formalism that we shall use for our lattice treatment in section 4 



One may gauge the above group 577(2) and arrive at the continuum action ( |3.44| ), which 
as we discussed above describes the low-energy continuum field theory limit of the large U 
Hubbard model (|2.29| ),( p73~5"l) . In this way, as we shall discuss below, one can generate the 
fermion condensate A3 dynamically. In this context, energetics prohibits the generation 
of a parity-violating gauge invariant SU(2) term 
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and so a parity-conserving mass 
term necessarily breaks |L2] the SU(2) group down to a T3 — U(l) sector 0, generated by 
the (T3 Pauli matrix in two-component notation. 



3.3 Connection with Superconductivity 

We now compare the model presented in this article and that of ref . || , which is known 
to exhibit unconventional parity invariant superconductivity, upon coupling the system to 
external electromagnetic potentials A^. First we note that there is an important physical 
difference between the two models, concerning the mechanism for mass generation. In our 
model in this article the gauge group that generates dynamically the fermion mass term is 
the strongly-coupled statistical Us(l), while the r 3 — £/(l)-remnant of the weakly-coupled 
SU(2) group is weakly coupled, and as such incapable of inducing mass generation. On the 
other hand, in ref. |J the fermion gap that lead to superconductivity was due to the t% — 
U(l) gauge boson. This may lead to important differences between the finite-temperature 
phase-diagrams of the two models. Such studies are left for future investigations. 



Nevertheless, as far as the mechanism of superconductivity is concerned, the two models 
appear to be qualitatively similar, and it is in this sense that the large-spin treatment 
of ref. pi is justified by the results of the present work. Indeed, the global U em (l) 



symmetry, which is a subgroup of the local symmetry of the ansatz (|2.17|) , corresponds 
to the electromagnetic symmetry in the statistical model. This symmetry can be gauged 
by coupling the action (|3.43|) to an external electromagnetic field on the spatial plane as 
in ref. [HI. 



5 On the continuum, of course, \£ = ^^70, with 70 a 2 x 2 Dirac matrix, and the hermiticity properties 
of the bilinears depend on the representation of the Clifford algerba chosen JlTJ] . 
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Figure 1: Anomalous one-loop Feynman matrix element, leading to a Kosterlitz-Thouless- 
like breaking of the electromagnetic U em (l) symmetry, and thus superconductivity, once 
a fermion mass gap opens up. The wavy line represents the SU(2) gauge boson B^, 
which remains massless, while the blob denotes an insertion of the fermion-number current 
J M = l I / 7 M \I / . Continuous lines represent fermions. 



As discussed there, then, superconductivity is obtained upon the opening of the gap in 
the fermion (hole) spectrum due to the one-loop anomalous effect corresponding to the 
following Feynman matrix element, depicted in fig. 1: 

S a =<B a \J u \0>, a = l,2,3; J„ = ¥7^ (3.54) 



with \l/ four-component spinors, which correspond to the continuum limit of ( 2.37|) . It 
should be stressed that as a result of the colour group structure only the massless B^ 
gauge boson of the SU(2) group, corresponding to the 03 generator in two-component 
notation, contributes to the graph. The result is [0, |J: 

S =< Bl\J v \0 >= (sgnM) W -^ (3.55) 



where M is the parity-conserving fermion mass (or the holon condensate in the context of 
the doped antiferromagnet). This observation is consistent with the symmetry-breaking 
patterns of the U em (l) group since the B^ colour component remains massless, and there- 
fore plays the role of the Goldstone boson ||. As discussed in ref. |7|, |[, this unconven- 
tional symmetry breaking however does not have a local order parameter, and thereby 
resembles, but is not identical to, the Kosterlitz-Thouless mode of symmetry breaking 



The massless Gauge Boson B^ of the unbroken U(l) subgroup of SU(2) is responsible for 
the appearance of a massless pole in the electric current- current correlator ||, which is the 
characteristic feature of any superconducting theory. In this sense, in ref. the field B^, 
or rather its dual <ft defined by d^cf) = e^pduB^, was identified with the Goldstone Boson 
of the broken U em (l) (electromagnetic) symmetry. In the non-Abelian context there are 
also Goldstone bosons associated with the breaking of the SU{2) symmetry [ 12| . These 
will be discussed in the next subsection. 
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4 Dynamical Gauge Symmetry Breaking on the Lat- 
tice 



In this section we derive the symmetry breaking patterns, and discuss, in detail, the 
excitation spectrum of the theory obtained from the effective long-wavelength lattice ac- 
tion ( |3.43| ). We are interested in the effective action of the holon degrees of freedom, after 
integrating out the fractional-statistics Us(l) field. From the above discussion it becomes 
obvious that this field plays an auxiliary role in the spin-separation ansatz, and as such 
it should be integrated out in the effective action of the physical degrees of freedom. 



We shall concentrate on the (3\ = strong coupling limit for the Us {I), which from the 
point of view of the doped Hubbard model corresponds to an infinite-?/ limit. In this limit 
the Us (I) gauge field may be easily integated out in the path integral with the result (T! 



J dVd^d^exp(-S eff ) (4.56) 



where 

Seff = fa £(! - trV p ) + £ lnJo(V^) 



y ifl = K^^Vi^i+^i+MV^t (4.57) 

and Iq is the zeroth order Bessel function. The quantity may be written in terms of 
the bilinears 

MabaB = ^i,b,f3^i,a,on a, 6 = colour, a, (3 = Dirac, i = lattice site (4.58) 
The result is: 

Vi » = -KV[M«(-7^M^)( Tm )^ (4.59) 

In the analogue language of particle physics |12[ the quantities Mw would represent 
physical meson states. In the context of our spin-charge separation ansatz the mesons 
would be composite states of holons. We have already seen that the physical electrons are 
composites of magnon-holons. In the theory (|3.43|) the magnon degrees of freedom have 
been integrated out. In this context, the low-energy (long- wavelength) effective action is 
written as a path-integral in terms of gauge field and meson states |l2fl 



Z = J[dVdM}exp(-S eff + £trlnM (i) ) (4.60) 

i 

where the meson-dependent term in ( fOKj ) comes from the Jacobian in passing from 



fermion integrals to meson ones p3 . 
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In ref. |12] a method was presented for identifying the symmetry-breaking patterns 
of the gauge theory ( |3.43| ), by studying the dynamically-generated mass spectrum. The 
method consists of first expanding J2i,n ^hiy/Ui^) in powers of y^, and concentrating on 
the lowest orders, which will yield the gauge boson masses, whilst higher orders describe 
interactions. Keeping only the linear term in the expansion yields |~2] 



m/ (v« - \vi, = -^ 2 tr[M«(- 7/ ,)^M^( 7/ ,)^] (4.61) 

It is evident that symmetry-breaking patterns for SU(2) will emerge out of a non zero 
VEV for the meson matrices AfW. 



Lattice simulations of the model ( |3.43| ), with only a global SU(2) symmetry, in the 



strong Us(l) coupling limit (3\ = 0, and in the quenched approximation for fermions, have 
shown [|35] that the states generated by the bilinears Ai and A2 (c.f. ( |3.51| )) are massless, 
and therefore correspond to Goldstone Bosons, while the state generated by the bilinear 
A3 is massive. In the context of our statistical model (c.f. fl2.39| )) these meson states may 
be expressed in terms of the holon operators as: 



Ai = -i(*i* 2 - *a*i)i = -2i(^i^2 - iptip 



1 i 



A 2 , = (W 2 + = -2(^2 + 1>Ui)i (4.62) 

and the bilinear A3 is given by (|2.40|) 

A 3 , = 1 - * 2 tf 2 )i = -2(ViVi - V4V> 2 ); (4.63) 

The fact that members of the triplet SU(2) representation acquire different masses is 
already evidence for symmetry breaking. We shall confirm this explicitly later on. For 
the moment we note that lattice analyses |JS| show that in the strong coupling limit 
Pi = the condensate u =< .A3 > and the mass of A3 are infinite. Of course the masses 
and the condensate are finite for finite fli, which is the case of finite-C/ Hubbard models 
(c.f. ( |2.34|) ). In addition, in this approximation this is the only meson state that develops 



a non-zero vacuum expectation value (VEV). This therefore constitutes a prediction for 
the infinite U Hubbard model and the spin-separation ansatz ( 2.17| ). The fact that the 



VEV of the Goldstone Boson states Al,2 vanish implies the absence of a 'spin-flip' (on 
the average) at a site: < ^|^ 2 ,i >=< 4'l,ii J i,i > = 0> which is also consistent with the 
slave-fermion constraints ( [2.2C| ). This is also comforting from the point of view of the 
equivalence of the above U — > 00 Hubbard model with that of ref. 0, whose symmetry 
breaking dynamical patterns are characterized by the absence of a local order parameter 0. 

6 The absence of VEVs for the Goldstone Bosons Ai,2 eliminates a potentialy dangerous source of 
a possible appearance of a local order parameter in the model. Notice that the dynamical breaking 
of the electromagnetic U em symmetry as a result of the holon condensate occurs without a local order 
parameter H. 
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One has the following expansion for the meson states in terms of the 577(2) bilinears 



M (i) = ^ 3 (0 ff 3 + Ax{i)a x + ^2(^)^2 + At(i)l + 

i[B^ + B hli {i)^a x + B 2 ^a 2 + £3,^7^3] (4-64) 

with fi = 0, 1, 2, 7^ are (antihermitean) Dirac (space-time) 2x2 matrices, and o"i, % = 1, 2, 3 
are the (hermitean) 2x2 SU(2)-'colour' Pauli matrices. Note that the VEV of the matrix 
< >= ua 3 is proportional to the chiral condensate. Upon substituting ( |4.64|) in 

(|4.61|) , taking into account that the SU(2) link variables may be expressed as: 

V ¥ = cos{\B ifl \) + ia.B iti sin(\B ifl \)/\B iti \ (4.65) 

and performing a naive perturbative expansion over the fields B one finds: 

\nl (y/y~) oc K 2 u 2 [(B}^) 2 + (B 2 ^) 2 ] + interaction terms (4.66) 

From this it follows that two of the SU(2) gauge bosons, namely the B l ,B 2 become 
massive, with masses proportional to the chiral condensate u: 

B ,2 boson masses cx K 2 u 2 (4.67) 

whilst the gauge boson B 3 remains massless. 

This mass term breaks 577(2) to a £7(1) subgroup, and in view of the above analysis one 
recovers the effective action for the massless modes occuring in the large-spin treatment 
of ref. ||, and reviewed in section 2. It is understood that a full analysis for finite values 
of j3\ is necessary, before definite conclusions are reached in connection with the exact 
properties and physical implications of the ansatz (|2.17| ) for finite U doped Hubbard, or 



t — j, models. We hope to come back to these issues in the future. 



We would like now to draw the reader's attention to the similarity of the above mech- 
anism for symmetry breaking with the situation in the adjoint gauge-Higgs model p7 . 



There, the SU(2) symmetry is also broken down to a £7(1) whenever the constant multi- 
plying the Higgs-gauge interaction is larger than a critical value. In our case the role of 
this constant is played by K 2 , as can be seen by the formal analogy between the adjoint- 
Higgs-gauge interaction terms and ( |4.61| ). Of course, in our approach symmetry breaking 
was achieved due to the infinitely strong Us(l) coupling. In view of the above analogy 
with the adjoint-Higgs model |37H , however, one may speculate that interesting phase di- 
agrams for the symmetry breaking of SU(2) could also emerge due to the K 2 coupling, 
in a way independent of the £7s(l) coupling. In this respect, we would like to stress once 
again that in the context of our statistical models the amplitude K is proportional to 
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the doping concentration in the sample, K oc Jrj. Since the adjoint-Higgs-like symmetry 
breaking requires strong enough coupling, then the above analysis, if true in this context, 
may be seen to suggest a natural and simple explanation -in the context of a gauge theory 
- of the fact that in planar antiferromagnetic models of finite- [/-Hubbard or t — j type, 
antiferromagnetic order is destroyed, in favour of superconductivity, above a critical dop- 
ing concentration. As mentioned at the end of section 2, this point of view seems to be 
supported by preliminary results of lattice simulations |yj . More detailed investigations 
along this line of thought are left for future work. 



5 Conclusions and Outlook 

In this article we have discussed lattice models for planar spin-| Heisenberg Antiferro- 
magnets away from half filling (doped). We have worked in the infinite U — > oo limit of 
the Hubbard model, which is characterized by the Gatzwyler projection, namely a con- 
straint of no more than one electron per lattice site. Upon implementing a spin-charge 



separation ansatz (|2.17|) , in a way consistent with holon spin flip, we have argued that 
the doped model is still characterized by a local SU(2) x Us(l) x U em (l) symmetry upon 
coupling to external electromagnetic fields. Of these, the Us(l) is an auxiliary 'statistical' 
gauge symmetry, associated with the fractional statistics of the spin and charge excita- 
tions in the ansatz ( |2.17| ). This possibility arises because of the planar spatial structure 
of the lattice model. 



We have argued that for strong enough Us{l) couplings, dynamical generation of a holon 
condensate can occur, with the result of breaking the SU(2) group to r 3 — U(l). This is the 
same local phase symmetry as the one characterising superconducting effective theories of 
doped antiferromagnets in large-spin S — > oo treatments || |J, although the mechanisms 
for mass generation are different. Nevertheless, the superconductivity scenaria appear 
qualitatively similar. In this way we have explained two things in a dynamical way: (i) 
the breaking- as a result of doping- of the local SU(2) spin symmetry that characterizes 
half-filled large-t/ Hubbard models, and (ii) the qualitative justification of large spin 
treatments and in particular the suppression of intrasublattice hopping of holes. Indeed, 
the latter is associated with massive SU(2) gauge boson states, which acquire their masses 
through holon condensation. There are many features of the models that still have to be 
worked out. Finite-f/ treatments and extension of these ideas to t — j models are worth 
pursuing. Given the dependence of the coupling constants of such models on the doping 
concentration in the sample, then a renormalizat ion-group study of the respective phase 
diagrams could provide useful quantitative information on the order of magnitude of the 
maximum doping concentration for superconductivity, and, in general, shed more light 
on the physics of the spin-charge separation in the models. We hope to come to a more 
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systematic study of such issues in the future. 



Further consistency checks of our approach may also come from a study of the renor- 
malization group structure of the normal phase of the model in the infrared. By normal 
phase we mean the phase where there is no dynamical opening of a gap. In this respect 
we mention that in three space-time dimensions the natural coupling constant appearing 
in the Lagrangian of a U(l) gauge theory with fermions is a parameter with dimensions of 
y / mass. In analytic Schwinger- Dyson treatments one can define a dimensionless coupling, 
which is essentially the ratio of the coupling constant over a characteristic mass scale of 



the theory, playing the role of the ultraviolet cut-off [gj • in a recent series of papers |38 
it was argued that this dimensionless coupling decreases slowly with the momentum scale. 
Its growth towards the infrared regime, however, is cut-off by the appearance of a non- 
trivial infrared fixed point. The latter phenomenon is responsible for deviations from 



fermi- liquid behaviour [^, |38f] , and - if the infrared fixed-point value of the coupling is 
strong enough [27|- also for mass generation. These features are expected to persist in 
the present model. However, in the present case, the full non-Abelian SU(2) x Us(l) 
symmetry will be present in the normal phase. A full analysis along the lines of ref. f38fl 
remains to be done. 



Above we have dealt with relativistic low-energy limits, obtained by linearizing about 
specific points on the fermi surface for the holons. As argued in ref. this may still 
capture certain qualitative features of realistic non-relativistic holon models. Eventually, 
one would like to be able to extend quantitatively the above results to non relativistic cases 
as well. We mention, however, that our relativistic limits may be related to condensed 
matter systems with fermi surfaces that have nodes. Such systems are known to exist 
in nature, and in particular they are antiferromagnetic planar systems with a strong 
spin-chain anisotropy as far as Heisenberg interactions are concerned f\. Upon doping 
and linearization around holon-fermi-surface nodes, one might then obtain the effective 
relativistic models discussed in this work and in ref. |J. 

An important issue we would like to raise as a result of the present work is the fact that 
non-Abelian local gauge symmetries, arising in the strong U Hubbard antiferromagnets, 
imply the possibility of existence of non perturbative effects (monopole-instantons in 
the form of Hedgehog configurations etc). Their precise role in the superconductivity 
mechanism associated with these models needs to be investigated in detail |], jl3| . This 
becomes particularly important in view of the claimed association of this scenario for 
superconductivity with Kosterlitz-Thouless-like phase transitions ||. There are important 
similarities between the two scenaria, since both are characterized by the absence of local 

7 We thank A. Tsvelik for useful information on the existence of such materials. 
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order parameters for the Goldstone bosons associated with the symmetry breaking. It 
is known that in Kosterlitz-Thouless transitions the symmetry breaking occurs when 
non-perturbative degrees of freedom are liberated. A preliminary analysis || [HJ in the 
effective theory model of ref . || , which, as a result of the present work, may be viewed as 
an effective theory of the massless degrees of freedom of the non-Abelian case, has shown 
that non-perturbative effects appear to be bound in pairs in the superconducting phase. 
This issue deserves however further investigations that require going beyond perturbation 
theory. 



In this latter respect, we mention that the treatment of non-perturbative effects requires 
exact results. Of course the superconductivity mechanism advocated in ref. || occurs 
through an anomaly, which is an exact one loop result. However, this is not sufficient for 
an exact quantitative treatment of the low-energy effective action. However, it is known 
that exact results in effective action treatments in higher than one spatial dimension can 
be derived in certain supersymmetric non-Abelian gauge theories, as a result of special 
non-renormalization theorems and strong- weak coupling duality symmetries ||40|| . In such 
theories, one invokes a duality symmetry to map a strongly-coupled problem to a weakly- 
coupled dual model which can be solved exactly. 



We now remark that t — j models are known, under certain restrictions among their 



parameters - namely t = j, to exhibit hidden super symmetries in space time [41]. There 



are graded algebras among the three possible states on a lattice site of the t—j model 
| a >= {|0 >,|1 >, 1 2 >}, corresponding to the empty, spin up and spin down states re- 
spectively. The model is supersymmetric up to a shift in the chemical potential, in the 
sense that there exist two supercharge operators , a = 1,2 (SU(2) 'spin' index), con- 
necting fermi and bose sectors and leaving the action invariant. So far this supersymmetry 
structure was not given any dynamical significance. This is because this supersymmetry 
refers to electron operators. Our ansatz ( |2.17|) , however, which implies electron substruc- 
ture, when and if extended to this case, might imply hidden supersymmetries among 
holon and spinons. These might have non-trivial consequences on the dynamics, follow- 
ing the spirit of ref. POJ, provided one could extend it to this case. In such a context, 



the superconductiviity model of ref. could be viewed as an effective theory of the light 
degrees of freedom, arising in the gauge-symmetry-broken phase of a supersymmetric 
SU(2) x U(l) x U em {l) field-theory model of a doped antiferromagnet with t = j. 



At present, we lack any microscopic dynamics underlying ( |2.17| ) that would allow us to 
check on its generalization to the t = j case and on the existence of the above-conjectured 
supersymmetric structure. At any rate, we believe that our work of associating holon 
condensation with a dynamical breaking of a Yang-Mills gauge theory in doped antifer- 
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romagnetic planar systems is an interesting observation, which deserves further serious 
investigations. We do hope to come back to a study of some of the above-mentioned issues 
in due course. 
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